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ABSTRACT 

The orbital properties of infalling satellite halos set the initial conditions which control 
the subsequent evolution of subhalos and the galaxies that they host, with implications 
for mass stripping, star formation quenching, and merging. Using a high-resolution, 
cosmological A^-body simulation, I examine the orbital parameters of satellite halos as 
they merge with larger host halos, focusing primarily on orbital circularity and peri- 
center. I explore in detail how these orbital parameters depend on mass and redshift. 
Satellite orbits become more radial and plunge deeper into their host halo at higher 
host halo mass, but they do not significantly depend on satellite halo mass. Addi- 
tionally, satellite orbits become more radial and plunge deeper into their host halos 
at higher redshift. I also examine satellite velocities, finding that most satellites infall 
with less specific angular momentum than the host halo virial value, but that satellites 
are 'hotter' than the host virial velocity. I discuss the implications of these results to 
the processes of galaxy formation and evolution, and I provide fitting formulas to the 
mass and redshift dependence of satellite orbital circularity and pericenter. 

Key words: methods: numerical - galaxies: halos - galaxies: kinematics and dynam- 
ics - cosmology: theory. 



1 INTRODUCTION 

In hierarchical structure formation, dark matter halos grow 
via accretion of both diffuse matter and virialized satellite 
halos. Beyond driving just halo mass growth, the nature of 
satellite accretion governs the subsequent evolution of both 
the satellite and host halos (hereon, 'satellite' and 'host' 
refer to the lower and higher mass halo, respectively, dur- 
ing a merger) . Since g alaxies for m at the centers of dark 
matter halos jWhite fc Roes 19?!; iBlumenthal et al.|[l98^ : 
iDubinskil [1994I : IMo et all i1998| '). the dynamics of infaUing 
satellite halos will also influence the evolution of the accom- 
panying galaxy populations. 

After infall, a satellite halo can survive as a substruc- 
ture halo (subhalo) of a larger host halo ( Tormen et al.l 
1 19981 : iKlvpin et ail Il999l : iMoore et all 1 19991 ). As it or- 
bits, a subhalo continues to host a (satellite) galaxy un- 
til it tidally disr upts or merges with the host halo's cen- 
tral galax^Q (e . g.,ISpringel et al.|[200ll : 'Kravtsov et al."2004l: 



Satellite orbits therefore influence the nature of 
galaxy mergers and the galaxy merger rate. Galaxy 
mergers are expected to drive morphological evolution 
(Toomre & Toomre 1972; Hausman fc Ostrikei)[l97i : IWhitel 
^1_978; Barnes fc H ernquist 1996,), and the properties of the 
post-merger gala:xy depend o n the relative dynamics of 
the galaxies du r ing the merger (jBoylan-Kolchin et al.l l2005l : 
ICox et al]|2006l : iHopkins et al.|[2008l ). It is thus necessarv to 
have cosmological predictions for how close satellite galaxies 
come to the central galaxy during first pericentric passage 
and the relative velocities distributions when they merge. 

After infall, satellites also can experience mass loss 
from tidal stri pping , tidal heating, and d isk shocking (e.j 



g 1 

Ostriker et al] Il972l: [Cnedin et al.' 1999*; 'P ekel et alj|2003l : 



Tavlor fc Babulll2004l : lb'Onghia et al..,201Q ). Satellite galax- 



ies are also thought to experience truncated star forma- 
tion and/or morphological evolution, arising from ram- 
pressure stripping, adiabatic heating (stran gulation), ha- 



Zentner et al.l l2005l: IConrov et al] [2OO6; Nataraian et at. 
20091 : IWetzel fc White! I2OIOI ). A satellite galaxy's survival 



rassment , and / or tidal shock heating (e.g., Gunn fc Gott 



timescale depends on its orbital parameters at the time of 
accretion, with galaxies on highly circular orbits surviving 
longer than those that rapidly plunge toward halo center. 



^ Additionally, ^ 10% of sate llite galaxies merge with another 
satelli te within the host halo jAngulo et al. I I2OO9I : I Wetzel et al.l 
l2009bl lah. 



1 197^ : iMoore et al.ll 19981 : lAbadi et al.H 19991 : iMcCarthv et al 
200§). Resonant stripping processes may also drive the 
evolution of dwarf s pheroidals around the Milky Way 
dP'Onghia et al.ll2009l ). The efliciency of all of these pro- 
cesses depends critically on the details of satellite galaxy 
orbits. 

Many semi-analytic models of satellite galaxy evolution 
contain prescriptions fo r dynamical friction survival times 
and tidal stripping (see iBaughl 12006] . for a recent review). 
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A number of authors have provided detailed fits to satelhte 
surv ival times that depend on a satelhte's circula rity at in- 
fall l|Bovlan-Kolchin et al.ir2008l : |jiang et al.ll2008l '). All these 
models require proper initial conditions of satellite orbits as 
a function of both halo mass and redshift to be an accurate 
depiction of galaxy evolution in a cosmological context. 

Several authors have examined the orbital param- 
eters of infalling satellite halosin simulati ons llXormenl 
19971: IVitvitska et al.l[200^: lBenson|[20 05: Wan g et al.ll2005l : 



Zentner et al.l l2005l : Khochfar fc BurkertI 12009 ') However, 
most work has focused on orbits only at z ~ 0, and lim- 
ited statistics and dynamic ranges have inhibited a robust 
investigation into possible mass and redshift dependence. 
The nature of large-scale structure formation suggests that 
satellite accretion will depend on halo mass and redshift. Tri- 
axial collapse models of halo formation predict that more 
massive halos form from more spherical regions, with ac- 
creting matter containing less specific angular momentum 
JZerPovichl Il970l: iBardeen et al'l 1 19861 : lEisenstein fc Loebl 



Sheth et al 



1 12001]) 



This is reflected in the nature of 
which low-mass halos reside primar- 



1 19951 : 

the cosmic web, in 
ily along/within filaments while hig h-mass halos resid e at 
the intersection of several filaments (|Bond et al.lll993 ). In 
this picture, matt er infall onto mass i ve halos occurs via 
narrow filaments jColberg et all Il999l : lAubert et all |2004 
iFaltenbacher et aLlbOOsF T and is more radial than matter in- 
fall onto low-mass halos, which are comparable in size to 
their filament and thus experience more isotropic infall. 

The nature of mass accretion may also vary with red- 
shift, with implications for the formation of the earliest 
galaxies. For instance, hydrodynamic simulations suggest 
that galaxy formation at z > 2 proceeds through highly 
radia l flows of co ld gas that penetrate deep into the host 
halo (|Keres et al ]l2005; Dckcl et al. 2009), behavior not ob- 
served in the local Universe. While these radial flows are 
influenced by complex gas processes, the dynamics of mat- 
ter at infall are governed primarily by gravity, so this trend 
suggests that mass accretion onto massive halos at high red- 
shift is significantly more radial than at low redshift. 

In this paper, we use a high-resolution A'^-body simula- 
tion of cosmological volume to examine the orbits of satellite 
halos at the time of infall into larger host halos. In addition 
to exploring the distributions of satellite orbital parameters, 
we also examine how the nature of satelhte accretion varies 
with halo mass and redshift. The combination of high res- 
olution and large volume allows us to study the nature of 
accretion from dwarf galaxy masses (10^° M©) to mas- 
sive galaxy clusters (10^^ M©) with good statistics from 
z = to 5. 



2 METHODS 

2.1 Simulations & Halo Tracking 

To find and track halos, we employed a dissipationl ess A ^- 
body simulation using the TreePM code of .Whitei (|2002D . 
This simulation used ACDM cosmology (Jim = 0.25, Qa = 
0.75, h = 0.72, n = 0.97 and erg = 0.8) in agreement 
with a wide a rray o f observations JPunklev et al.l 120091: 



Kowalski et"all l2008l : IVikhlinin et all l2009l : iPercival et al.l 
20ld ). For high mass and force resolution, the simulation 



evolved 1500'^ particles in a 200/i~^Mpc box, with a par- 
ticle mass of 1.64 x 10* Mq and a Plummer equiva- 
lent smoothing of 3ft~^kpc. Initial conditions were gener- 
ated using second-order Lagrangian Perturbation Theory at 
z = 250 where the RMS was 20% of the mean inter-particle 
spacing. 45 outputs were stored evenly in ln(a) from z; = 10 
to 0, with an output time spacing of ~ 6 50 Myr at z = 0. 
Note that the same simulation was used in IWetzel fc Whit3 
(|2010h . 

Halo s were found using the Friends-of- Friends (FoF) al- 
gorithm ijPavis et al.l Il985l ) with a linking length of b = 
0.168 times the mean inter-particle spacinglf] We kept all 
FoF groups with more than 60 particles, corresponding to 
10^'' ^"^Mq. This resolution level ensures that FoF halo 
masses are accurate to within 10% () Warren et al.l liood l. 
Merger trees were constructed from the set of halo catalogs 
by specifying a parent-child relationship. To be a 'parent', a 
halo must contribute more than half of its mass to a 'child' 
halo at the next simulation output. Thus, a parent halo can 
never have more than one child. A merger was then identi- 
fied as a child halo with more than one parent. Note that 
only distinct FoF halos (no subhalos) were used in this work. 



2.2 Ejected Halos &: Re-mergers 

Once a satellite halo has fallen into a larger halo, it can 
retain its identity as a bound subhalo. In some cases a sub- 
halo can become ejected from its host halo, either from 
being on an initially unbound orbit or from a scattering 
event within the host halo. Up to 50% of all satellite halos 
within 2 — 3-Rvir of a host halo are a recently ejected popu- 
lation, particularly at low satellite halo mass (since subha- 
los experience severe mass stripping as they pass through 
a host halo) and for satellites with highly eccentric o rbits 
(|Gin et al.ll2005l : iLudlow et al.ll2009l : fWang et al.ll2009l ). In- 
cluding these re-merging satellite halos would bias the or- 
bital distribution results both by double counting a single 
satellite halo across redshift and by artificially enhancing 
the population of highly elliptical orbits. 

Since our merger tree does now allow a parent halo to 
have more than one child, recently ejected halos are nec- 
essarily parent-less. Thus, we identified ejected halos (and 
their subsequent re-mergers) by finding halos which have no 
parent halo but had at least one of their 20 most bound par- 
ticles within some halo at the previous output We define 
infall as the first time a halo merges with a more massive 
halo, and we used this event to compute satellite orbits, 
discarding subsequent re-mergers. A re-merger was not dis- 
carded, however, if the satellite halo had more than doubled 
its mass since becoming ejected, since in this case we con- 
sidered a new halo to have formed. 



2.3 Calculating Orbits 

To calculate the orbital parameters of satellite halos at infall 
into a more massive host halo, we used the simulation halo 



^ The longer linking length of 6 = 0.2 is often used, but it is more 
susceptible to joining together distinct, unbound structures. 
^ We checked that this method g ives consistent re s ults a s com- 
pared with the subhalo catalogs of lWetzel fc Whitei l|2010l '). 
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merger trees to find halos which were about to merge, that 
is, become joined within an isodensity contour by the FoF 
algorithm. We then computed the orbital parameters using 
halo masses, positions, and velocities in the output prior to 
merging. Positions and velocities were those of the halo's 
most bound particle, which is expected to correspond to 
a hosted galaxy. Velocities were calculated using physical 
coordinates, including Hubble flow. 

With these halo properties, we calculated orbital pa- 
rameters by treating the two halos as isolated point parti- 
cles in the reduced mass frame (limitations of this approx- 
imation are discussed below and in the Appendix). Given 
satellite mass Msat and host halo mass Mhost, with reduced 
mass ^ — A/sat Mhost / (Afsat + A/host), satellite separation 
r = Thost — fsat and velocity r = f-host - rsat, the orbital 
energy is 

GAfsatA^host 



P 1 -2 



and the angular momentum is 



(1) 



L = iiv XV . (2) 
With these, we computed the orbital eccentricity. 



e = 4 1 + 



(GAfsatA4oet)V 



pericentric distance (from host halo center), 



(1 -I- e)GA4"sat A/host M 



and apocentric distance 



(3) 



(4) 



(5) 



(1 - e)GAfsatAfhcst/i 

Orbital circularity is defined as the ratio of the orbit angular 
momentum to that of the circular orbit with the same energy 
and can also be related to eccentricity 

m 



ME) 



(6) 



In the point particle two-body approximation, the 
above quantities are all conserved throughout the orbit, so 
they do not depend on the separation of the satellite from 
the host halo. Two of the above quantities (for example, cir- 
cularity plus pericenter) are sufficient to uniquely describe 
an orbit. An alternate description can be given by the satel- 
lite's radial, Vi, plus tangential, Ve, velocity, though they 
need to be defined at a given radius, which here will be at 
the host halo virial radius. To do this, we used the conser- 
vation of orbital energy and angular momentum to evolve 
the velocities from the satellite's measured location to where 
its center crosses the host halo virial radius, i?vir, which was 
derived from the h alo FoF mass and co ncentration assuming 
a spherical NFW (|Navarro et al.|[l996l ) density profile. 

Hereon, all satellite halo distances and velocities are 
scaled to the host halo virial radius and virial circular ve- 
locity, Kir = Ay/GA^virZ-Rvir. 



2.4 Calculating Orbital Distributions 

The aim of this work is to examine the orbital parameter 
distributions of satellite halos as they first pass through the 



virial radius of a larger host halo. However, because of the 
finite time spacing of simulation outputs, satellite halos were 
identified at a variety of distances from the host halo virial 
radius prior to merging. This presents some difficulty in di- 
rectly measuring orbital parameter distributions, since satel- 
lite halos can experience mass loss (or growth), dynamical 
friction, and tidal forces, all of which alter the orbital pa- 
rameters as satellites traverse significant distances. 

To circumvent this issue of finite time resolution, we 
estimated satellite orbital parameter distributions by select- 
ing only satellite halos that were within a small separation 
of the host halo virial radius at the output prior to merging 
(as d one by IVitvitska et al.ll2002l : iBensonllioOsI : IWang et all 
l2005h . Specifically, we found satellite halos that were about 
to merge with a larger host halo and whose edge (given by 
its own virial radius) was in range [1, 1 + Ar] R^ir of the 
virial radius of the host haloQ We used Ar = 0.25, finding 
no significant difference using smaller values. 

Since more radial orbits spend less time within the given 
radial shell, to properly estimate the overall orbital distri- 
butions one must scale the orbital counts by the crossing 
time within the shell. Specifically, we scaled the counts by 
Af/fcross, where At is the output time interval and tcross is 
the orbital crossing time for a satellite halo edge to orbit 
from a separation of Ar to 00 

Since this method of estimating satellite orbital distri- 
butions considers only satellites close to the location of in- 
terest, it is little-affected by mass stripping, dynamical fric- 
tion, and tidal forces, and the error on the estim ated or- 
bital parameters is expected to be less than 10% (|Bensonl 
l2005h . There are, however, limitations to the point particle 
two-body approximation, which ignores multi-body interac- 
tions and extended, triaxial halo profiles (see the Appendix). 
Some satellites (up to 20% at 2 = 0) had their rperi or rapo 
within the radial selection shell, so these objects did not have 
a well-defined crossing time even though they did merge with 
the host halo in the next output. This population was dom- 
inated by satellites on highly circular orbits, such that rapo 
was never beyond the outer radial shell. In weighting their 
counts for orbital distributions, these satellites were given a 
shell crossing time equal to the longest well-defined crossing 
time satellite in the given output, which amounts to a mini- 
mal weighting without discarding them. Additionally, some 
merging satellites were on formally unbound orbits, though 
this represents less than 2% of satellites regardless of mass 
or redshift0 



* We used the location of the satellite edge since the FoF algo- 
rithm merges halos when the satellite edge crosses the host halo 

5 This 



corr ection for crossi n g tim e was employed by 

BensonI ||2005|) and IWang et al.l | |2005| '| but was neglected in 



Vitvitska et al.l (120021 ) 



° Another method to estimate satellite orbital parameter distri- 
butions is to consider all satellite halos which are about to merge, 
regardless of their distanc e from the host halo in the o utput prior 
to merging, as done by iKhochfar fc Burkerj 1 I2OO6I ). However, 
they found that in order to select systems which conserve en- 
ergy and angular momentum, they were restricted to 'isolated' 
merger systems which do not change in total mass by more than 
10% between outputs. Indeed, we found that using all satellites 
which are about to merge yields as high as 15% of orbits being 
formally unbound, since satellite halos on highly eccentric orbits. 
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Figure 1. Distributions of circularity, r], and pericenter, rperi 
(scaled to host halo virial radius), of infalling satellite halos at 
z = 0, for halos > 10^'^ Mq. Errors bars indicate Poisson error 
in each bin. Vertical lines show average circularity and median 
pericenter. Dotted curves show fits to the distributions (Eqs. [7] 
andig. 



Finally, we also examined how using Af200crit instead of 
FoF mass influences the resultsQ We found that the averages 
of the orbital distributions from the two methods are within 
error, and thus the specific halo finding algorithm used does 
not sig nificantly affect the results, in agreement with similar 
tests of lBensoni l|2005l '). 
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Figure 2. Distributions of tangential, Vg, radial, Vr, and total, 
Vtot , velocity (scaled to host halo virial velocity) of infalling satel- 
lite halos at 2 = 0, for halos > 10^" hr^ Mq. Vertical lines show 
average ve locities. Dashe d an d dot-dashed curves show distribu- 
tions from lBensonI l|2005h and lWang et al.l l|2005h . 



3 ORBITAL DISTRIBUTIONS AT Z = 

We begin by examining the distributions of satellite orbital 
parameters for all infalling satellites at z = 0. Figure[T]shows 
the distribution of circularity and pericenter, scaled to the 
host halo virial radius, for all halos > 10^'' /i^^ Mq. Cir- 
cularity shows a broad distribution, peaking at 77 = 0.52 
(which corresponds to an ecce ntricity of e = 0.85), in good 
agreement with earlier work (|Tormen 19971: IZentner et al.l 
I2OO5I : IWang et al.ll2005l : iKhochfar fc Burkertll2006l ). Satel- 
lite orbits tend to be neither highly radial nor highly cir- 
cular. Pericenter shows a distribution falling exponentially 
with radius, with a median value of rpcri = 0.21 and falling 
to zero for rperi > 1. Most infalling satellites are on orbits 
taking them close to the center of the host halo. 

While an orbit can be fully classified by its circularity 



which are more likely to be found at large radii in the output 
prior to merging, become significantly more bound prior crossing 
the host halo virial radius. We found that while using this alter- 
nate method yielded quantitatively different orbital distributions, 
it yielded qualitatively similar results regarding the orbital dis- 
tribution shapes and their mass and redshift dependence. 
^ The relation of FoF(b = 0.168) mass to spherical overden- 
sity (200crit) mass depends on halo concentration and redshift, 
though for the regimes considered here they are within ~ 15%. 



and pericenter, another such combination is radial and tan- 
gential velocity as measured at the host halo virial radius. 
Figure [5] (top and middle) shows these satellite velocity dis- 
tributions, scaled to the host halo virial velocity. Since the 
velocities are computed at the host halo virial radius, the 
tangential velocity is equivalent to the satellite specific an- 
gular momentum: L/Lvir = Vs/Vvir. Most satellites infall 
with less specific angular momentum than the virial value 
of the host halo, with an average Ve = 0.64. However, this 
does not mean that satellite mergers do not contribute sig- 
nificantly to the angular momentum build-up of halos, since 
the a ctual spin of a halo is only ~ 5% of its virial value 
(e.g., IVitvitska et~al] I2OO2I ) . Indeed, the merging of satel- 
lites is the primarily driver of a halo's angular momentum 
(|Vitvitska et al.ll2002l ). 

As opposed to satellite tangential velocity, radial veloc- 
ity is typically comparable to the host halo virial velocity, 
with an average K = 0.89. Taken together, these imply that 
most satellites are infalling 'hotter' than their host halo. 
Figure [2] (bottom) demonstrates this explicitly, showing the 
satellite total velocity distribution, again scaled to the host 
halo virial velocity. On average, satellites are infalling with 
~ 15% higher velocity than the matter within the host halo. 

To compare with previous work. Fig. [2] als o shows 
f its to the velocity distribu tions at z = of iBensonI 
(2005") and lWang et"all (|2005h . which represent a Gaussian 
distribution for radial velocity and approximately a two- 
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Figure 3. Distributions of circularity and pericenter at 2 = 0, for 
satellite halos of mass 

mass ranges. Vertical lines show average circularity and median 
pericenter. Dotted curves show fits to the distributions (Eqs. [7] 
andlSjl. Satellite orbital distribution shapes depend on host halo 
mass. 



dimensional Maxwell-Boltzmann distribution for tangential 
velocity. While the radial velocity distribution seen here ex- 
hibits a somewhat narrower profile, the distributions show 
overall broad agreement, particularly in the averages and 
peak locations. Some of the differences may be attributed to 
previous wo rks using somewhat different cosmology and that 
IWang et al.l |2003) selected satellite subhalos (halos which 
may have already merged) within a radial shell. More im- 
portantly, the previous works examined satellite orbits at 
a higher mass regime, and as will be explored in the next 
section, satellite orbital distributions depend on halo mass. 

Finally, the above tangential velocity distribution ex- 
plains why pericenter exhibits an exponential distribution. 
At fixed mass ratio, rpori oc Vg /{I + e) (Eq. |4|. Since eccen- 
tricity is essentially limited to < e < 1, the rperi distri- 
bution is dominated by the Vg term. Vg is described by a 
two-dimensional Maxwell-Boltzmann distribution, P{Vg) ~ 
Vee"'^""^"-' , so under transformation, the distribution of 
Vg, and hence rpcri, is exponential. Note that considering 
orbits in extended host halo profiles will modify this expo- 
nential distribution somewhat (see the Appendix). 



4 MASS DEPENDENCE 

We next explore how satellite orbital distributions depend 
on halo mass at 2 = 0. We first examined the dependence on 
host halo mass by selecting satellite halos in the mass range 
-[^gio.0-10.5 1^-1 -^^^ used the large dynamic range of 

the simulation to explore the dependence across four decades 
in halo mass. 

Figure [3] shows the distribution of circularity and peri- 
center, as in Fig. [T] but for two host halo mass ranges. The 
distributions clearly change shape, being skewed to both 
lower circularity and lower pericenter for more massive host 



I I I I llj 1 1 — I I I I ll| 1 1 — I I I 1 1 1| 1 1 — I I I 1 1 L 
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Figure 4. Dependence of satellite average circularity and me- 
dian pericenter on host halo mass at 2 = 0, for satellite halos 
of mass iqIO O^iO-S Mq. Dotted curves show standard errors 
of the mean/median. Inset panels show same, but dashed curves 
indicate the standard deviation (or 68 percentile width) of the 
full distribution. Satellite orbits are both more radial and plunge 
deeper at higher host halo mass. 



halos, which drives the average/median of the distributions 
(vertical lines) down. 

Figure |4] demonstrates explicitly the dependence of av- 
erage circularity and median pericenter on host halo mass. 
Circularity shows no dependence up to ~ 3 x 10^^ Mq, 
but above this mass satellite orbits become less circular with 
increasing host halo mass, with average circularity dropping 
nearly 20% across the mass range. Interestingly, the turnover 
corresponds to the value of Af, , the characteristic halo mass 
scale of collapse, at z; = 00 Median pericenter decreases 
more strongly with host halo mass, falling by more than a 
factor of 2 across the mass range and showing no rollover at 
low mass. Overall, satellite orbits are both more radial and 
plunge deeper into their host halo at higher host halo mass. 

Similarly, Fig. [S] (top and middle) shows the host halo 
mass dependence of satellite average tangential and ra- 
dial velocity. Both velocity components remain below the 
host halo virial velocity at all mass scales. Tangential ve- 
locity monotonically declines with host halo mass, falling 
by 30%, implying that satellite accretion contributes angu- 
lar momentum less efficiently onto higher mass halo. Ra- 
dial velocity declines rapidly with host halo mass up to 
^ 3 X 10^^ Mq, beyond which it remains nearly flat. 

These trends explain the host halo mass dependence of 
circularity and pericenter. Below ~ 3 x 10^^/i^^Mq, both 
Vg and K decline with host halo mass, giving rise to con- 
stant circularity. At fixed circularity (eccentricity) and Msat , 
rpcri/-Rvir OC Vg / M^^^^ (Eq. O, SO rpori coutiuues to decline 



* More specifically, Mt{z) is the mass at which a{M,z), the 
variance of the linear power spectrum at redshift z smoothed on 
scale M, equals the threshold for linear density collapse, Sc = 
1.69. 
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Figure 5. Dependence of satellite average tangential, radial, and 
total velocity on host halo mass at z = 0, for satellite halos of 
mass 10l° °-l°-5 /i-l M0. Inset panels show same, but dashed 
curves indicate the standard deviation of the full distribution. 
The degree to which satellite infall velocities are biased relative 
to their host halos depends on host halo mass. 



rapidly with mass. Above ~ 3 X 10^^ h'^ Mq, declining Ve 
and constant K cause both circularity and pericenter to de- 
cline with mass. 

Additionally, Fig. [5] (bottom) shows the host halo mass 
dependence of satellite average total velocity. While satel- 
lite infall is usually 'hotter' than the host halo virial veloc- 
ity, the magnitude of this effect decreases with host halo 
mass, falling by 20% over the mass range. Interestingly, this 
leads to a crossover such that satellites infalling onto halos 
> 10^* Mq are instead on average 'colder' than the host 
halo. This mass trend is driven by the interplay of halo mass 
with its environment. High-mass halos generically dominate 
the local potential field, and so satellites naturally fall in 
with velocity comparable to the halo virial velocity. By con- 
trast, many low- mass halos reside in proximity to higher 
mass halos which significantly boost the local potential field, 
causing satellites to fall into low- mass halos faster th an their 
virial velocity (see also Fig. 8 in lWetzel et al. I l2008l . on how 
mass, separation, and environment affect halo mergers). 

To highlight the significance of these mass dependencies 
with respect to the overall orbital distributions, the inset 
panels of Figs. [4] and [5] show the same as the main panels, 
but the dashed curves indicate the full distribution width 
(standard deviation). While mass dependence is significant 
in each case, note that the mean/median changes by at most 
one standard deviation across the mass range probed here. 
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Figure 6. Dependence of satellite average circularity and me- 
dian pericenter on satellite halo mass at z = 0, for host halos of 
mass 10^*~^^ Mq. No significant, systematic dependence on 
satellite mass is found at any mass scale or redshift. 



While the above results demonstrate that satellite or- 
bital parameters depend on host halo mass. Fig. [6] explores 
whether circularity and pericenter depend on satellite halo 
mass, for host halos of fixed mass IO^-^'^/i-^Mq. While 
Fig. [6] suggests a possible drop in circularity and pericenter 
for very massive satellites, the large scatter across satellite 
halo mass does not lead to any clear, systematic trends. 
Moreover, we found no clear trends with satellite mass for 
any orbital parameters, including varying host halo mass 
or redshift. Thus, we conclude that the nature of satellite 
infall is controlled by host halo mass and is little affected 
by satellite mass. This supports the physical picture that 
orbital dynamics are governed by the most massive halo 
within a region, and that less massive satellites effectively 
act as massless tracers of the potential field. 



5 REDSHIFT EVOLUTION 

We now turn to explore the dependence of satellite orbital 
parameters on redshift. Since the halo mass function de- 
clines significantly with redshift, examining the mergers of 
all halos above the resolution limit (> 10^° Mq) would 
convolve any redshift dependence together with the mass 
dependence of the previous section. So to isolate redshift 
trends, we examined halos in a fixed mass range: host halos 
of mass 10^2.0-12.5 ^-i (corresponding to M* halos at 
2 = 0) and satellite halos of mass 10^°"" Mq. 

Figure [7] shows the distribution of circularity and peri- 
center at three redshifts. Similar to the dependence on host 
halo mass in Fig. [3l the distributions are skewed to both 
lower circularity and lower pericenter at higher redshift. Fig- 
ure [S] demonstrates more explicitly the dependence of aver- 
age circularity and median pericenter on redshift. Average 
circularity falls nearly 30% from z = to 5, while pericenter 
falls more rapidly to less than half of its 2: = value. Over- 
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Figure 7. Distributions of circularity and pericenter at different 
redshifts, for host lialos of mass 

10^2.0-12.5 1^-1 a^j^jj satellite 
halos of mass 10^^~^^ Mq. Dotted curves show fits to the 
distributions (Eqs. [7]and|8ll. Satellite orbital distribution shapes 
depend on redshift, though the distribution widths remain broad 
at all redshifts. 
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Figure 8. Dependence of satellite average circularity and median 
pericenter on redshift, for host halos of mass 1012-0—12.5 ^-1 jyj^ 
and satellite halos of mass lO^^^^i /i^^Mq. Satellite orbits are 
both more radial and plunge deeper at higher redshift. 



all, satellite orbits become more radial and plunge deeper 
into their host halos with increasing redshift. 

Figure [5] also shows the redshift evolution of satellite 
velocities. Tangential velocity declines with redshift, imply- 
ing that satellite accretion less efficiently contributes angu- 
lar momentum growth at higher redshift. Conversely, radial 
velocity increases with redshift, approaching the host halo 
virial value at z ~ 5. Interestingly, the contrasting evolu- 
tionary trends of radial and tangential velocity lead to no 
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Figure 9. Dependence of satellite average tangential, ra- 
dial, and total velocity on redshift, for host halos of mass 
10^2.0-12.5 j^-l a^jjj satellite halos of mass lO^O"" h'^ Mq. 
Tangential velocity declines with redshift, while radial velocity 
increases, leading to nearly constant total velocity. 



significant evolution in satellite total velocity, which at all 
redshifts remains 'hotter' than the host halo virial velocity. 

Note that the widths of the total distributions do not 
evolve significantly with redshift, retaining standard devia- 
tions of ~ 0.2 as shown in Fig. [T] Thus, as with mass de- 
pendence, the redshift evolution of the mean/median out to 
2; = 5 is comparable to the distribution width itself. 

What drives the redshift evolution of satellite orbits? 
One possibility is that redshift dependence is simply a man- 
ifestation of the trends with mass from the previous sec- 
tion. Analytical triaxial collapse models of halo for mation 
(JJardccn ct al. 1986; Eisciistcin & Locb 1995; Shc th et al.l 
I2OO1I ) predict a self-similarity in the nature of matter in- 
fall with redshift at fixed M/Mt{z), since it is at M« that 
halos transition from being located along filaments to be- 
ing at the intersection of several filaments. In this picture, 
the above redshift dependence is driven by probing higher 
Ahost /M*{z) at higher z. 

To investigate this, we selected host halos of mass 
Mhost /Mt{z) = 270 at each redshift (corresponding to 
10" ''^ h-^ Mq at 2 = down to 10^"'^ h'^ Mq at z = 3) 
and satellite halos in the range 10^" °"^°-^ /i"^ M©. Fig- 
ure [To] (top) shows that satellite infall into host halos of 
constant Ali^ost /M,{z) indeed exhibits nearly constant cir- 
cularity. However, pericenter (bottom) is not constant but 
instead increases with redshift by a factor of 2 (opposite to 
the trend at fixed mass). Furthermore, at fixed Mhost /A'/* (z) 



© 2010 RAS, MNRAS OOO.mfTTI 



8 Wetzel 



0.55 




Figure 10. Same as Fig.[8l but for host hales of mass 270 M, (z) 
(10" '^5 /i-l M0 at 2 = down to 10^°-^h-^MQ at z = 3) 
and sateUite halos of mass iqIO O-IO.S f^-l jy^^ When scaled by 
M*(z), circularity remains nearly constant while pericenter in- 
creases significantly. 




Figure 11. Distributions of circularity and pericenter at differ- 
ent redshifts, for host halos of mass 270 M, (z) and satellite halos 
of mass lol'' 0~l'^-3 ft-l Mq. Dotted curves show fits to the dis- 
tributions (Eqs. [3 and |8]l. The circularity distribution is nearly 
universal for host halos of fixed Afijost/A^* (z)- 



we found a significant increase in radial, tangential, and to- 
tal velocity with redshift. Thus, while the redshift evolution 
can be partially understood simply as a manifestation of 
mass dependence, intrinsic redshift dependence does exist. 







Pi 


7i 




Co 


3.38 


0.567 


0.152 
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Ci 


0.242 


2.36 


0.108 
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Ro 


3.14 


0.152 


0.410 


(l + z)-4 


Ri 


0.450 


-0.395 


0.109 


(H-z)-4 



Table 1 . Fit parameters for d, Ri, given by Eq. [g] for use 
circularity (Eq. [tJ and pericenter (Eq. [S]! distributions. 



6 FITS TO ORBITAL DISTRIBUTIONS 

We now seek to provide analytical fits to the orbital dis- 
tributions of the two parameters of primary interest, circu- 
larity and pericenter. The results of the previous sections 
show that the shapes of these orbital distributions depend 
sensitively on both mass and redshift. However, Fig. 1101 
shows that, at fixed Mhost/Mt{z), circularity remains nearly 
constant with redshift. Figure [11] (top) demonstrates this 
more explicitly, showing that the overall circularity distribu- 
tion remains approximately universal at fixed Mhost/Af*(z). 
However, the distribution of pericenter is neither redshift- 
invariant at fixed Mhost (Fig. [7} nor at fixed Mhost/Af*(z) 
(Fig. EI. 

The orbital distributions of satellite circularity, rj, and 
pericenter, rpcri, are given to good approximation across 
mass (10^""^'^ h-^ Mq) and redshift (2 = - 5) by the fol- 
lowing simple functional forms 



=Co(Mhost,z)r?'°'(l 

ari 



n) 



Ci(Mi„rt,z) 



(7) 



d/ 



dro 



i?o(Mhost,2)exp{-[rperi/i?i(MhoBt,2)]°-'*^} . (8) 



Parameters Ci (Mhost, z) and i?i (Mhost, z) describe the mass 
and redshift dependence and take a similar functional form 

li 



Mh, 



M^z) 



(9) 



where the values for ai,Pi,^i and the function gi{z) are 
given in Table [T] Since the circularity distribution is nearly 
constant when scaled by Mhost /M«(2), gi[z) — 1 for d. 
While the pericenter distribution does not exhibit such red- 
shift invariance, the scaling (7i(2)Mhost/Mt(2) with gi{z) = 
(1 + 2)"* for Ri fully encapsulates the redshift dependence. 

Since parameters Co (Mhost, 2) and i?o (Mhost, «) are 
merely normalizations for probability distributions, their 
mass and redshift dependencies are implicitly given by those 
of Ci (Mhost, 2) and i?i (Mhost, 2). We include fits to their 
mass and redshift dependencies for completeness. 

Finally, given the assumed cosmology the evolution of 
M, (2) is approximated to within 5% up to 2 = 7 by 



log [M,(z)/h~'^ Mq] = 12.42 - 1.562 + 0.0382^ 



(10) 



The orbital distribution fits given by Eqs. [7] and [S] pro- 
vide a good description of both the mass and redshift de- 
pendence of satellite orbits. The accuracy of these fits is 
demonstrated explicitly in Figs. [T] O [71 and llll 

We note two caveats to these orbital parameter fits. 
First, these are valid for satellites at the time of infall and 
were derived assuming point particle orbits. Thus, they are 
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likely to be modified well after infall as satellites experi- 
ence extended, triaxial halo potentials (see the Appendix), 
as well as dynamical friction, mass stripping, and scatter- 
ing processes. Second, while we did not examine in detail 
correlations between the circularity and pericenter distri- 
butions, some correlation is expected, with more circular 
orbits tending t o have larger pericenters, though with sig- 
nificant scatter llTormenlll997l : iGill et al.ll20oi : lBenson|[2005l : 
iKhochfar fc Burkertll2006l '). We leave a more detailed inves- 
tigation into these behaviors to future work. 



7 SUMMARY & DISCUSSION 

We used a high-resolution A'^-body simulation of cosmo- 
logical volume to track the mergers of halos from dwarf 
galaxy masses (10^° M©) to massive galaxy clusters 
(10^^ M-q) across a large redshift range [z — Q to 5). 
We explored the orbital parameters of satellite halos at in- 
fall, when they cross within the virial radius of a larger host 
halo. The main results are as follows: 

• Satellite orbital parameters exhibit broad distributions, 
with typical standard deviations of ~ 0.2. For all resolved 
halos at 2 = 0, the average values of the distributions are: 
circularity t) = 0.52, (median) pericenter rperi/J?vir = 0.21, 
tangential velocity Ve/V^ir = 0.64, radial velocity V^/Vvir = 
0.89, and total velocity Kot/Vvir = 1.15. 

• Satellite orbits are more radial and plunge deeper into 
their host halos at higher host halo mass, but orbits are not 
significantly affected by satellite halo mass. 

• Infalling satellites are typically 'hotter' than the host 
halo virial velocity, except onto the most massive host halos 
where satellite are slightly 'colder'. 

• At fixed halo mass, satellite orbits become more radial 
and plunge deeper into their host halos at higher redshift. 

• The satellite circularity distribution exhibits almost no 
redshift evolution for host halos of fixed Mhost/Af* (z), im- 
plying that redshift dependence of circularity at fixed mass is 
simply derivative of varying M/Mt(z). However, pericenter 
exhibits significant increase with redshift at fixed A/host /Af*. 

As explored in ^ the orbital distributions we find at 
z = when stacking halos of all masses agree well with 
those of previous work. However, previous results of possi- 
ble ma ss and/o r reds hift dependence of satellite orbits are 
mixed. iBensonI (|2005l ) found evidence that satellite orbits 
become more radial at higher halo mass scales but was un - 
able to quantify this f urther, while IVitvitska et al.l (|2002l ) 
and lWang et"aLl (|2005l ) saw no such halo m ass dependence 
over a limited mass range. By contrast, IVitvitska et al.l 
(2002") found that satellite angular momentum decreases 
with increas i ng sa tellit e mass to host halo mass ra tio, while 
IWang et all (120051 1 and lKhochfar fc BurkertI ||2006D found no 
dependence on mass ratio, though again oyer a li mited mass 
ratio range. Finally, the results of iBensonI ll2005ll su g gested 
trends with redshift while those of IVitvitska et all i|2002l l 
did not. In most cases, previous work was limited in terms 
of merger statistics (in some cases, examining infall into a 
handful of halos) and dynamical range (unable to explore 
both M > M. and M < M*). The results here demon- 
strate clear dependence of satellite orbits on host halo mass 
and redshift and no significant evidence for dependence on 



satellite mass. This is broadly consistent w ith the predic- 
tions of a nalytical triaxial collapse models (Ba rdeen et al.l 
Il98d : ilisenstein fc Loeb .1995, : .Sheth et al. 200li), in which 
more massive halos arise from a more spherical Lagrangian 
volume with less angular momentum. We emphasize that the 
existence of mass and redshift dependence implies that fits 
to orbital dist ributions based on stacking halos of all masses 
at z ~ fe.g.. lBenso nl l2005l : [Zentner et al. l l2005l : IWang" et al.l 
l2005h are not universally accurate. 

The mass and redshift trends seen here have implica- 
tions for various aspects of galaxy formation and evolu- 
tion. For example, recent work suggests that galaxy for- 
mation a t z > 2 proceeds through narrow streams of 
cold gas (jKerei et al.ll2005l : IPekel et all 120091 ') . fundamen- 
tally different behavior than seen in the local Universe. 
iDekel fc Birnboiml (|2006l ) argued that this difference arises 
in part because massive (~ 10^^ h"^ Mq) halos at high red- 
shift form at the intersection of narrow filaments, while at 
low redshift such halos are more likely embedded within a fil- 
ament and experience wide-angle inflow. Our results on the 
redshift evolution of satellite accretion qualitatively support 
this picture, but it is not clear that the evolution seen here 
is strong enough to imply a fundamental change the nature 
of accretion a.t z ^ 2.5, when most orbital parameters dif- 
fer from their z — values by ~ 20%. This suggests that 
the above results are driven more strongly by evolving gas 
physics than by the nature of mass accretion. 

Our results also have clear implications for satellite 
galaxy evolution within groups/clusters. For example, us- 
ing the parametrizati o n of satellite survival time given in 
iBovlan-Kolchin et al.l (|2008l ). changes in circularity across 
mass and redshift from the global z = value can lead 
to reductions of satellite lifetimes of 30% or more. Addi- 
tionally, environmental effects such as ram-pressure strip- 
ping of gas and tidal stripping of stars a re expected to oc- 



cur primarily at orbi t al pericenter (e.g., Dekel et al.l |2003| : 



iTavlor fc Babulll2004l : iMcCarthv et al.ll2008t ). Beyond their 
dependence on evolving gas physics, our results suggest that 
satellite galaxy quenching and morphological evolution pro- 
ceed more efficiently and rapidly at higher group/cluster 
masses and higher redshift. 

While this work focuses on the orbital parameters of 
satellite halos at the time of infall, it is not immediately 
clear how well these orbital distributions and their mass and 
redshift dependencies persist to satellite populations well af- 
ter infall, as the orbits become affected by dynamical fric- 
tion, tidal stripping, and extended halo potentials (see the 
Appendix for some estimate of the latter). There has been 
some work examining satellite subhalo orbits within hos t 
halos (|Gill et al.1 12001 : iReed et~alll2005l : [Sales et al.ll2007l ). 
which found orbital parameter distributions similar to those 
here (for example, rj ~ 0.5). This lack of circularity evolution 
after infall is plausible since dynamical friction is expected to 
be i nefficient in altering orbital parameters su ch as circular- 
ity (I van den Bosch et al.|[l999l ). Furthermore. iFaltenbaclieil 
(201^) recently found that the orbits of satellite subhalos 
within host halos at 2: = are more radially biased in more 
massive host halos, which suggests that the results here re- 
main valid well after infall. We will pursue a more robust 
investigation into the evolution of satellite orbits after infall 
in future work. 
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There is also possi ble observational evidence in support 
of the trends seen here. iHerbert-Fort et"al] (|2008| ) examined 
the orbital velocities of galaxies in local galaxy clusters, find- 
ing highly asymmetric velocity distributions consistent with 
satellites large ly retaining their infalling orbits. Promisingly, 
iBiviano &: Po ggianti ( 2009) examined satellite galaxy orbits 
in galaxy clusters from z — to z — 0.8 and found evidence 
that satellite orbits are indeed less isotropic (more radial) 
at higher redshift. 

Finally, the results on satellite velocities also have impli- 
cations for relating satellite dynamics to those of the overall 
group/cluster. At lower host halo masses, satellite velocities 
become significantly 'hotter' than than the host halo, imply- 
ing a possible systematic biasin g in using satellite velocity 
dispersions to infer halo masses (|Faltenbacheij|2010l . found a 
similar trend with mass for the velocity bias of satellite sub- 
halos within host halos). Future work will involve a more 
detailed analysis of satellite velocity bias and its mass and 
redshift dependence. 
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Figure Al. Distribution of pericenter (scaled to host halo 
virial radius) for infalling satellites halos at 2 = 0, for all halos 
> IO^^/i^^Mq. Solid curves shows distribution derived by nu- 
merically integrating satellite orbits after infall assuming a spher- 
ical NFW density profile for the host halo. For comparison, dotted 
curve shows distribution derived via the point particle approxi- 
mation, as in Fig. [T] Vertical lines show corresponding median 
values. 

over at small radii. Overall, the spherical NFW pericenter 
distribution is skewed to higher values, causing the median 
to increase from 0.21 (point particle) to 0.32 (NFW). 

Thus, the pericenter distributions given in this paper 
may moderately underestimate the true distances of closest 
approach that satellites experience. However, beyond this 
systematic offset, we find that using an NFW profile pro- 
duces no significant change in the mass and redshift trends 
seen above. Furthermore, other factors such as triaxial pro- 
files, dynamical friction, mass stripping, and scattering pro- 
cesses will also alter a satellite's orbital properties after in- 
fall, and future work will explore these in more detail. 



APPENDIX A: IMPACT OF HALO PROFILE 
ON ORBITS AFTER INFALL 

Satellite orbital circularity and pericenter as explored in this 
work were derived from the point particle approximation 
and so are only strictly valid at the time of infall. Circu- 
larity and pericenter are not generally conserved for orbits 
within more realistic extended halo potentials. In particu- 
lar, the actual pericentric distance that a satellite experi- 
ences after infall can differ systematically from the value 
derived at infall. Here, we briefly explore the impact of an 
extended host halo profile on the satellite pericenter distri- 
bution at z = 0. To do so, we proceeded in the same manner 
as above (outlined in , except that we modeled host halos 
as having spherical NFW profiles as given by their measured 
concentrations, and we numerically integrated the orbits of 
(point particle) satellites to obtain their distance of closest 
approach. 

Figure IATI shows this pericenter distribution at z = 0, 
for all halos > 10^" M©. For comparison, the dotted 
curve shows the original distribution that was obtained from 
the point particle approximation (same as in Fig.[l]). Because 
an extended halo profile contains less mass at smaller radii, 
a satellite experiences reduced gravitational focusing as it 
nears halo center, causing the pericenter distribution to roll 
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